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The field of quantum computation currently lacks a formal proof of experimental feasibility. 
Qubits are fragile and sophisticated quantum error correction is required to achieve reliable quantum 
computation. The surface code is a promising quantum error correction code, requiring only a 
physically reasonable 2-D lattice of qubits with nearest neighbor interactions. However, existing 
proofs that reliable quantum computation is possible using this code assume the ability to measure 
four-body operators and, despite making this difficult to realize assumption, require that the error 
rate of these operator measurements is less than 10 -9 , an unphysically low target. High error 
rates have been proved tolerable only when assuming tunable interactions of strength and error rate 
independent of distance, which is also unphysical. In this work, given a 2-D lattice of qubits with only 
nearest neighbor two-qubit gates, and single-qubit measurement, initialization, and unitary gates, 
all of which have error rate p, we prove that arbitrarily reliable quantum computation is possible 
provided p < 7.4 x 10 -4 , a target that many experiments have already achieved. This closes a 
long-standing open problem, formally proving the experimental feasibility of quantum computation 
under physically reasonable assumptions. 



A 2-D array of qubits with tunable nearest neighbor in- 
teractions is a believable experimental target [IHI]. The 
surface code [5HTU] can be implemented optimally using 
such an array. Qubits undergoing very general unwonted 
evolution can be accurately modeled as suffering random 
Pauli X and Z errors [IT] . This is a reasonable model as 
quantum error correction (QEC) involves repeated oper- 
ator measurements that project noisy qubits onto states 
that differ from the desired state by Pauli operators. An 
expandable pattern of operator measurements M{ asso- 
ciated with the surface code is shown in Fig. [TJ The 
operators Mj are called stabilizers [12) . 

Let the state of the qubits in Fig.[l]be |\&) and, with- 




FIG. 1: A small surface code. Larger codes can be constructed 
by expanding the pattern. Circles represent qubits. Each 
Mi represents an operator (tensor product of Pauli X or Z 
operators) that is measured to detect errors. Note that all 
operators commute. 



out loss of generality, let us assume that we have mea- 
sured all stabilizers Mj and observed the +1 eigenstate, 
so Mi |*) = |*) for all i. Note that [M^Mj] = for 
all i, j, so such a state exists. Suppose the central qubit 
suffers some general error E, where 

Eq. [2] shows that E is a linear superposition of no error 
(I), or an X, Z or XZ error. After the error E, a sub- 
sequent round of Mj measurements will project E to / if 
Mi = +1 for i = 3, 5, 6, 8, X if M t = +1 for i = 5, 6 and 
Mi = -1 for i = 3, 8, Y if M, = -1 for i = 3, 5, 6, 8, and 
Z if Mi = +1 for i = 3, 8 and M 4 = -1 for i = 5,6. All 
measurements Mj, i ^ 3,5,6,8, will remain +1. 

The earliest proof that arbitrary reliability could be 
achieved using the surface code assumed that each mea- 
surement Mi had a probability q of reporting the wrong 
result and that between each round of measurement each 
qubit suffered a Z error with probability p |13j . The proof 
showed that arbitrary reliability was achievable given suf- 
ficient qubits for p + q < 2.4 x 10 -11 . More recently, 
assuming perfect stabilizer measurements and a 3-D ar- 
ray of qubits each suffering a single error with probabil- 
ity p, assumptions that are equivalent to a 2-D array of 
qubits with periodic faulty measurements and qubit er- 
rors, arbitrary reliability was proved possible given suffi- 
cient qubits for p < 1.4 x 10 -9 [T3]. These error bounds 
are unphysically low. 

Error rates slightly above 10 -3 have been proved tol- 
erable in other quantum error correction codes when as- 
suming tunable interactions between arbitrary pairs of 
qubits with interaction time and error rate independent 
of the physical separation of the qubits [T51[TE] ■ However, 
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no physical machine possesses such interactions. As such, 
there is currently no formal proof that arbitrarily reliable 
quantum computation is experimentally feasible, raising 
serious questions about the validity of quantum compu- 
tation. In this work, we provide the much-needed formal 
proof of experimental feasibility. 

It has long been believed that error rates of 10 -4 are a 
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reasonable experimental target [T7l [18] . Experimentally, 
single-qubit measurement with error rate 10~ 4 and ini- 
tialization with error 10 -5 have been demonstrated [T§] . 
Single-qubit unitary gates have been demonstrated with 
error 2 x 10 -5 [2U]. Two-qubit gates are the most techni- 
cally challenging, with Bell state preparation with error 
7 x 10~ 3 the current state-of-the-art [21]. There is no 
physical reason to believe the technical challenges can- 
not be overcome and similarly low-error two-qubit gates 
achieved. All of these experiments were performed using 
ion traps [22], a technology well-suited to implementing 
a 2-D array of qubits with tunable nearest neighbor in- 
teractions [2j. We therefore seek a formal proof that ar- 
bitrarily reliable quantum computation can be achieved 
given a 2-D array of qubits with only nearest neighbor 
two-qubit gates, single-qubit measurement, initialization, 
and unitary gates, with all gates having an error rate 
p~ 10" 4 . 

The only nontrivial unitary gates we will use will be 
Hadamard H (H |0) = |+), H\l) = |-}) and controlled- 
NOT C'x 10) = |11), C x |ll) = |10)). Quantum 
circuits measuring the stabilizers Mi are shown in Fig. [2] 
Note that these circuits can be implemented in parallel 
across an arbitrarily large surface. A single round of 
error detection is defined to be a single parallel execution 
of these circuits. Initialization to |0) results in |1) with 
probability p, measurement reports the wrong eigenstate 
with probability p, H and identity I are followed by an 
error randomly chosen from X, Y, Z with probability p, 
and Cx is followed by an error randomly chosen from 
IX, IY, IZ, XI, ZZ with probability p. This error 
model is justified by the above discussion of arbitrary 
errors being projected to Pauli errors. 

Errors followed by even perfect Cx gates can result in 
multiple errors. For example, Cx |00) = |00), however if 
there is an X x error then C X X X |00) = C x |10> = |11) = 
ATiA^lOO). In words, the Cx copies X errors on the 
control qubit to the target qubit. Similarly, Cx |++) = 
|++), however if there is a Zi error then CxZ2 |++) = 

C x |H — > = I > = z i z 2 1++)- In words, the C x gate 

copies Z errors on the target qubit to the control qubit. 

Define a detection event to be a pair of sequential sta- 
bilizer measurements that differ in value. Fig. [3^, shows 
a single X error propagating through rounds of surface 
code error detection circuitry until two detection events 
are generated. X errors are detected by sequential pairs 
of Z stabilizer measurements, Z errors by sequential pairs 
of X stabilizer measurements. Given the right pattern 
of errors, any pair of sequential stabilizer measurements 
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FIG. 2: Quantum circuits measuring a. ZZZZ and b. XXXX 
operators. The |0) represents initialization, the H represents 
Hadamard, the Mz represents measurement of the operator 
Z, and the dot and target symbols connected by lines each 
represent a Cx gate. The Cx inverts the value of the tar- 
get qubit if the control (dot) qubit is in the state |1) and 
does nothing otherwise. For example, Cx(ot |0) + /3 |1)) |0) = 
a |00)+/3 jll). The interaction sequence is North, West, East, 
South, as shown in Fig. [3] 




FIG. 3: (Color online) 2-D surface code (grey). Time runs ver- 
tically. Squares represent initialization to |0) , circles represent 
Z basis measurement. Slashed squares represent initialization 
to |+), slashed circles represent X basis measurement (both 
achieved using Hadamard gates), a. A single error leading 
to a pair of detection events (long vertical ellipses). A detec- 
tion event is a sequential pair of measurements with differing 
value. Lines with arrows show the branching paths of error 
propagation, b. An error leading to a single detection event 
due to proximity to a boundary of the lattice. 



can be associated with a detection event. We associate a 
specific space-time location with each potential detection 
event, namely the space-time midpoint between sequen- 
tial measurements. 

For visualization purposes, we draw a sphere at every 
space-time location a detection event can occur. We draw 
a cylinder between every pair of space-time locations that 
can be associated with detection events generated by a 
single error. We call this structure of spheres and cylin- 
ders a lattice of dots and lines. Lattices are constructed 
by studying the propagation of all errors through the pe- 
riodic surface code quantum circuit V2'S \ 1241 . An example 
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FIG. 4: (Color online) Distance 4 example of a lattice of dots 
and lines with stochastically generated detection events. Dots 
(small spheres) correspond to space-time locations where the 
endpoints of error chains could potentially be detected. De- 
tection events (large spheres) correspond to space-time loca- 
tions where error chain end points have been detected. Lines 
link pairs of dots where a pair of detection events could poten- 
tially be generated by a single error. Darker lines link spatial 
boundaries to a single dot where a single detection event could 
be generated by a single error. 

of a surface code lattice is shown in Fig. |1J including de- 
tection events stochastically generated by errors during 
a simulation of the surface code. Note the lines leading 
to nowhere, indicating the nearby presence of potential 
errors leading to single detection events. This only oc- 
curs near boundaries of the array of qubits. Note also 
the regular yet nontrivial structure of the lines. 

There are two lattices, one associated with detecting X 
errors and one with detecting Z errors. The distance of a 
surface code is defined to the length in lines of the short- 
est path between disjoint boundaries of a lattice. Paths 
of lines connecting disjoint boundaries can be associated 
with chains of Pauli operators that commute with all 
stabilizers yet are not products of stabilizers themselves. 
Such chains of operators are called logical operators and 
manipulate the data stored in the surface code. 

For the purposes of this proof, we associate a weight of 
1 with each line. Given a lattice with random detection 
events, a minimum weight perfect matching is a set of 
paths through the lattice such that every detection event 
is incident on exactly one path, paths are allowed to ter- 
minate at boundaries, and the total weight of all paths 
is minimal. Edmonds' minimum weight perfect match- 
ing algorithm [25-25] can efficiently find such a set of 



paths. Other algorithms have been used to correct errors 
in the surface code [29ll32| . however matching currently 
remains the only algorithm that has been used to effi- 
ciently handle errors in implementations of the surface 
code making use of only two-qubit interactions and sin- 
gle qubit measurements. 

Define the line probability e to be the total proba- 
bility of all single errors leading to detection events at 
the endpoints of the line. Given Pauli error models for 
each quantum gate, with the strength of each error model 
scaled by a global parameter p, the probability of a given 
line can be expressed as a polynomial in p, with a linear 
relationship at low p [23] . 

Let Vbeannxnxn dot volume of lattice with non- 
cyclic boundaries. Assume each pair of opposing faces 
is a boundary of a distinct type. Let p be a character- 
istic physical gate error rate low enough to ensure the 
probability of any given line in V being associated with 
an error is less than some small e. For minimum weight 
perfect matching to fail to correct the errors in V, there 
must be at least one path of to > n + 1 lines connecting 
opposing boundaries containing at least [to/2] lines as- 
sociated with errors. A path corresponding to a logical 
error not containing at least [to/2] lines associated with 
errors can be matched with lower weight, contradicting 
the assumption of a minimum weight perfect matching. 

The number of paths of length to can be upper 
bounded by choosing three of the six faces of V as path 
starting points. This results in 3n 2 starting points. Each 
dot is connected to at most 12 neighboring dots (Fig. [2]). 
A general path from dot to dot that does not backtrack 
on itself must link to one of 11 of the surrounding 12 dots 
(the 12th is already part of the path). After the starting 
point is chosen, there are only to — 1 direction decisions 
to make. There are thus no more than 3n 2 ll m_1 paths 
of length to. Since we are only interested in non-self- 
intersecting, opposing boundary connecting paths, this 
is a generous upper bound. 

Given a particular path of length to, the probability of 
at least [to/2] of its lines being associated with errors is 




= 2 m e r ^ (7) 



The probability of a logical error is therefore no more 
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Provided e < 1/484, the above can be made arbitrar- 
ily small by increasing n. Since e can be expressed as a 
polynomial in p independent of n, this proves that using 
minimum weight perfect matching to correct errors in the 
surface code results in a finite threshold error rate. As- 
suming the surface code circuits and error models in |23j , 
in which the most errorprone line satisfied e < 14p/5, a 
lower bound to the threshold error rate of 7.4 x 10~ 4 is 
obtained. Given experimental achievements to date, this 
is sufficiently high to formally prove the experimental 
feasibility of arbitrarily reliable quantum computation. 

Eq. [12] also proves that logical errors are exponentially 
suppressed with code distance, implying extremely low 
logical error rates can be achieved with modest qubit 
overhead. Furthermore, as our proof is fundamentally 
based on minimum weight perfect matching and this al- 
gorithm is highly efficient [371 HE] > we have proved that 
the classical computing overhead is also modest. 
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